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Abstract 

Starting from approximate Skyrmion solutions obtained using the rational map 
ansatz, improved approximate Skyrmions are constructed using scaling arguments. 
Although the energy improvement is small, the change of shape clarifies whether the 
true Skyrmions are more oblate or prolate. 
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1 Skyrmions and Rational Maps 

The Skyrme model [1] is an effective field theory used to describe nuclei. For static fields, 
the classical solutions of the Skyrme model are stationary points (either minima or saddle 
points) of the energy functional 

e= i^L H tr _ ^ tr Rj]2) + m2tr(i ~ u) } d ^' (1) 

where x = (2:1,2:2, £3) axe the Cartesian coordinates, and Z7(x) is an S"?7(2)-valued scalar 
field, known as the Skyrme field, which describes pions nonlinearly. The gradient of the 
Skyrme field is captured by the current components Ri = dill U~ l (i = 1,2,3). In (JTJ), 
the energy and length units have been scaled away, leaving only the pion mass parameter 
m, which is the physical pion mass in scaled units. [In many past studies the pion mass 
term m 2 tr(I — U) was omitted.] In this study, when we consider massive pions, we assume 
that m=l, since only for m at or close to this value can the proton and delta masses be 
reproduced (for more details, see [2] and References therein). Also, when m = 1, nuclei are 
represented by Skyrmions that have a fairly uniform density, rather than a hollow, roughly 
spherical structure [31 H]. 

Finiteness of the energy requires that U approaches the identity matrix I at spatial 
infinity (and this boundary condition is also imposed if m = 0). Hence, M 3 is topologically 
compactified to a large 3-sphere. U is then a mapping from S 3 h-> SU(2), and is classified 
by the integer-valued degree (winding number) 

fl = --L/" tr(i?ii2 2J R 3 Kx, (2) 

which is a topological invariant. In particular, B classifies the solitonic sectors of the model 
and in applications to nuclei, B is identified with the baryon number of the Skyrme field 
configuration. A Skyrmion is a static soliton of minimal energy for given B, although 
sometimes, more loosely, a local minimum or saddle point of the energy functional is also 
called a Skyrmion. 

In [5], minimal energy Skyrmions with massless pions were approximated by an ansatz 
involving a rational map R between Riemann spheres. Here, a point in M 3 is labelled by 
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coordinates (r, z) where r is the radial distance from the origin and z = tan(0/2) e llf specifies 
the direction from the origin (9,<p are the usual spherical polar coordinates). Let R be a 
degree N rational map of the form R(z) = where p and q are polynomials in z with no 
common factors, such that max [deg(p), deg(g)] = N. Then the ansatz for the Skyrme field 
is 



where f(r) is a real profile function satisfying the boundary conditions /(0) = kn (for some 
integer k) and /(oo) = 0. It is straightforward to verify that the baryon number of this field is 
B = kN. In the remainder of this paper we only consider k — 1, thus B = N. An attractive 
feature of the rational map ansatz (j3J) is that it leads to a simple energy expression which is 
minimized by first minimizing with respect to the parameters in the rational map R(z), and 
then determining the profile function f(r) by solving an ordinary differential equation. This 
approach gives close approximations to true Skyrmions. Clearly, this procedure encounters 
difficulties when there are two or more Skyrmion solutions, either saddle points or genuine 
minima, with very similar energies. The ansatz might lead to the wrong energy ordering. 

For massless pions the minimal energy Skyrmions up to B = 22 are all well-approximated 
by the rational map ansatz [6] . However, the inclusion of the pion mass term in the energy 
functional makes the Skyrme model more realistic in describing real nuclei. The solutions 
change dramatically when m — 1, especially for large baryon number (B > 8). Some 
Skyrmions are known to be prolate (cigarlike) or oblate (pancakelike), for example, those 
with baryon numbers B = 8 and B = 12, respectively [I]. Spherical shell-like configurations 
given by the rational map ansatz are still useful starting points in the search for true solutions, 
but they are unstable [3j, because their hollow interiors have substantial potential energy 
when m — 1, and the fields relax to more compact structures, some being rather flat. This 
way the Skyrmion reduces its volume to surface area ratio and thus reduces its energy. 
Moreover, the Skyrmions with massive pions are smaller in size and exponentially localized, 
compared to Skyrmions with massless pions which are algebraically localized. This can be 
seen in Figure HJ 




(3) 



2 



2 Independent Length Rescalings 

Recently, Manton [7j studied the effect of independent length rescalings in the three Cartesian 
directions in the Skyrme model. The simple Derrick scaling identity [8] and further novel 
identities were derived, relating contributions to the total energy of a Skyrmion. For an 
exact Skyrmion solution, the identities are satisfied exactly. 

In what follows we will study the effect of scaling manipulations on the energy functional 
for various Skyrme field configurations which are not exact solutions, for both massless and 
massive pions. In particular, for an approximate Skyrmion described by the rational map 
ansatz, we can lower the energy by a rescaling of the form 

xi -> A1X1 , x 2 -> A 2 x 2 , £3 -> A3X3 , (4) 

with the parameters Aj not all the same. For a general Skyrme field, U(xi,X2, £3) is replaced 
under such a rescaling by U(xi, x 2 , £3) = U (XiXx, A 2 x 2 , X3X3) and the energy E of the rescaled 
field U(x) is the modified version of the original energy (JTJ), 

e - Mi- ^ - ^ w - - l¥* 

-i^itr ([R 2 , R s f) - iMltr ([R 3 , Rl f) + -j£-tr (I - U) )<f>x , (5) 

o A\ o A 2 AiA 2 A3 J 

where the current components Ri are evaluated for the original field U(pc). 

For an exact Skyrmion solution U(x), E is stationary with respect to A» at A» = 1 (since 
the energy is stationary with respect to any smooth change of the field that preserves the 
boundary condition). As in [7], we may restrict to special rescalings which together span all 
possibilities. In particular, a rescaling in the (2^2,^3) plane with no rescaling of x\ can be 
achieved by setting A 2 = A3 = A and Ai = 1 in (j3j). Then the derivative of E with respect 
to A vanishes at A = 1. This gives identity ([6]). Similarly, by permutation, the following 
identities are obtained: 
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We refer to the left hand sides of these identities as the scaling integrals {Xj : i — 1, 2, 3}. The 
sum of the above three identities is Derrick's identity for Skyrmions, which is also obtained 
by considering the uniform rescaling Ai = A2 = A3 = A. 

In [6], it has been shown that for m = and B > 7, there are increasingly many 
Skyrmions, with different symmetries and different shapes, whose energies are very close 
to the minimal value. The additional configurations are local minima or saddle points of 
the Skyrme energy, whose energies are difficult to distinguish numerically. In particular for 
B = 9 there are two Skyrmion possibilities with and symmetries; the first is probably 
the global minimum and the second a saddle point. For B = 10 there are at least four 
solutions close to minimal energy, but for B = 11 the minimal energy Skyrmion appears 
isolated. It can be difficult to ensure numerically that the fields have fully relaxed to a 
solution. One good way to test if they have would be to check whether the scaling identities 
(]6])-(IE]) are satisfied. 

The rational map ansatz gives good approximations to exact Skyrmion solutions, which 
generally do not satisfy all these scaling identities. An optimal rescaling should deform 
the rational map configurations closer to the exact Skyrmions. Moreover, after rescaling, 
the scaling identities will all be satisfied. The argument is as follows. We start with the 
approximate solution and evaluate the seven contributions to the energy that occur in ([Q), 
and with modified coefficients in fl5]). We then find the parameters Aj that minimize fl5]). 
This tells us how to rescale the initial field so as to obtain an improved approximate solution. 
The rescaled field satisfies the identities flS])-®, as it has minimal energy with respect to 
further rescalings. An alternative formulation of this, more convenient computationally, is 
to say that after rescaling, the identities 

X i: / \-l^tr(Rl) + l^tr([R 2 ,R 3 f) + ^tr(I-U)}d^ = 0, (9) 
J R 3 I I A 2 A 3 8 Ai AiA 2 A 3 J 

X 2 : \ \-\ A-tr (Rl) + I ^itr ([ife, Rtf) + (I - U) ) rf 3 x = , (10) 

J R 3 I Z A3A1 8 A 2 AiA 2 A 3 J 

[ \-l^tr(Rl) + l^tr([R 1 ,R 2 r) + -^tr(I-U)}d^ = (11) 
Jrz I 1 ^1^2 o A 3 AiA 2 A 3 J 

are satisfied, where Ai,A2,A3 are the parameters we find, and U, Ri, R 2 and R3 are the 

original field and current components before rescaling. 
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The energy is lowered by this rescaling, so one gets closer to an exact solution. More 
interesting, perhaps, than the small reduction of energy is the small change in shape. We 
learn from the calculation whether the true Skyrmion is more prolate or oblate than the 
approximate Skyrmion (which itself is rather round), or triaxial. 

Note that in all cases, the product of the rescaling parameters is very close to unity, i.e. 
A1A2A3 = 1, if the starting point is the optimised rational map ansatz with the profile function 
worked out numerically. This is because the starting point satisfies Derrick's identity. 

After rescaling a Skyrme field, the baryon density becomes 



where the right hand side denotes A1A2A3 times the original baryon density evaluated at 
x = (A1X1, A 2 x 2 , A3X3). Its integral is of course unchanged. In particular, if AiA 2 A 3 = 1, 
the baryon density transforms as a scalar quantity, so that the new density at x is the old 
density at x. 

In order to make some practical use of the scaling identities (E])-© it would be best to 
work with Skyrmions that are known to be, or are expected to be, very far from spherically 
symmetric. Examples are the Skyrmions with m = 1 and baryon number a multiple of 
four, composed of B = 4 subunits [I] as in the a-particle model of nuclei. These Skyrmions 
sometimes look like part of the infinite Skyrme crystal. 

We will calculate here the optimal rescaling of approximate Skyrmions given by the 
rational map ansatz, for B = 6, 9, 10 and 11, and for both m = and m — 1. In some of 
these cases the exact solutions are known, so we do not learn much new. For B = 9 and 
B = 11 the exact solutions for m — 1 are not known, so here we get some insight into the 
shapes of the true Skyrmions. Since the rescalings are all quite close to unity (within a few 
percent), it is important that the calculations of the parameters \ are not dominated by 
numerical errors in the energy integrals. In fact, the numerical errors in the rational map 
approach are smaller than 1%. 

We now describe our procedure for finding these configurations, and rescaling them. Ini- 
tially, a rational map is introduced with a given topological charge B and assumed symmetry. 
Then, the energy is minimized (using a simulated annealing process) with respect to the pa- 
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Figure 1: Baryon density isosurfaces of the rescaled (approximate) Skyrmions for B = 6, 9, 
10 and 11 Skyrmions with m = (first row) and m = 1 (second row). Each corresponds to 
a value of £> = 0.035. 
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rameters appearing in the rational map. In practice, this means minimizing an integral Z, 
given in [5], which depends only on the rational map R. Using this familiar approach, we 
have re-derived the results obtained in [5], [6] and [9]. The rational maps we obtain are the 
following: 

B = 6 

The rational map has symmetry, with the D\ generated by z i— > iz and z h-> 1/z, and 
has the functional form 

R(*)= 2 r + 4 + lV ( 13 ) 

z z [latz* + 1) 

The energy is minimized when a = 0.15853. The reflection symmetry arises because the 
parameter a is real. 

B = 9 

The energy minimizing rational map in this case has symmetry and has the form 

z (a + iPz* + z*) 
K[ > l + i(lz* + azs ' 1 ' 

where a = -3.37764 and (3 = 11.20311. 
B = 10 

The minimal energy Skyrmion for m = 1 is found to have D 2 h symmetry [5] . This symmetry 
is consistent with the a-particle model, where the field configuration consists schematically 
of a pair of cubic B = 4 Skyrmions separated by two B = 1 Skyrmions. The optimal rational 
map with this symmetry is of the form 

, a + /3z 2 + jz 4 + 5z 6 + ez 8 + z w 

{Z} ~ l + ez 2 + 5z A + 7 z 6 + (3z s + az 10 ' ( } 

with a = 0.2772, = -9.3594, 7 = 14.81, 5 = 4.977 and e = 3.015. 

B = 11 

The rational map here is Z^-symmetric, and of the form 

z 9 + az 6 + (3z 3 + ^ 
{Z> z 2 ( 7 z 9 + $z 6 + az* + 1) 1 ' 
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where a = -2.4719, = -0.8364 and 7 = -0.1264. 

For each of these optimised rational maps, the profile function f(r) is obtained by solving 

the radial equation 

,,, / 2Bsm 2 f\ 2f sin2f / . n N Xsin 2 f\ 2 
/" (l + ^T ± ) + "f + — F 1 (/ 2 - 1) - ""^J - sin / = , (17) 

using a shooting method. B is the baryon number and X is the integral mentioned above, 
evaluated on the optimised rational map. To avoid singularities at the origin, the so-called 
local analysis (see, for details, Reference [TO]) has been applied. Near r = 0, the profile 
function can be approximated by the Frobenius series f(r) = tt — ar a where the indicial 
exponent a is determined by taking the positive root of the quadratic equation a 2 + a — 2B = 
0, and a is the shooting parameter. The other boundary condition is /(oo) = 0. 

The Skyrme field is then transformed back from Riemann to spherical coordinates, i.e. 
(r, z, z) — > (r, 9, if), by inverting the equation z = tan (6/2) e llf , and the current components 
Ri are calculated. Numerical integrations, for evaluating the contributions to the energies E 
and E, and to the scaling integrals X, and Zj, are performed by the Gauss-Kronrod method. 

Finally the values of the rescaling parameters Aj are obtained by minimizing the energy 
(JSJ) using (once more) a simulated annealing process. They are presented in Tables [1] and |2j 
The values are slightly further from 1 for m = 1 than for m = 0. 



m = 


B = 6 


B = 9 


B = 10 


B = 11 


Ai 


0.98954 


1.01272 


1.02050 


1.00441 


A 2 


0.98954 


1.01272 


0.98811 


1.00441 


A 3 


1.02118 


0.97524 


0.99149 


0.99103 



Table 1: Values of the scaling parameters Aj for Skyrmions with massless pions. 

The energies of the approximate Skyrmions, before and after rescaling, are presented (up 
to five decimal places) in Table El while the values of the scaling integrals are presented 
in Tables H] and [5] and prove that after rescaling, the scaling identities are almost exactly 
satisfied. This is a check on our numerics. Note also that A1A2A3 is very close to 1, as 
anticipated. The changes in the energies are rather small, and much smaller than the 1% 
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m = 1 


B = 6 


B = 9 


B = 10 


B = 11 


\i 


0.98921 


1.01427 


1.02224 


1.00455 


A 2 


0.98921 


1.01427 


0.98785 


1.00455 


A 3 


1.02226 


0.97268 


0.99052 


0.99095 



Table 2: Values of the scaling parameters A, for Skyrmions with massive (m = 1) pions. 



B 


G 


E(m = 0)/B 


E(m = 0)/B 


E{m = 1)/B 


E(m = 1)/B 


6 


Du 


1.13726 


1.13675 


1.33151 


1.33092 


9 


Du 


1.11681 


1.11610 


1.31788 


1.31696 


10 


D 2h 


1.11218 


1.11172 


1.31561 


1.31504 


11 


D 3h 


1.10976 


1.10968 


1.31632 


1.31622 



Table 3: Values of energy per baryon before rescaling (TjQ) and after rescaling (J5]), for m = 
and m = 1. G is the symmetry group of the Skyrmion. 



777 — 


B = 6 


B = 6 


B = 9 


B = 9 


B = 10 


5 = 10 


B = 11 


B = 11 




0.00793 


-1.33 10- 7 


-0.00934 


1.59 10" 9 


0.00882 


1.29 10" 9 


-0.00332 


5.11 10" 8 




0.00793 


-1.05 10" 7 


-0.00934 


-3.90 10" 12 


-0.01509 


3.74 10" 9 


-0.00332 


-1.11 lO" 10 


2 3 


-0.01593 


1.21 10" 9 


0.01866 


1.04 10" 9 


0.00624 


8.44 10- 10 


0.00662 


-3.27 10" 11 



Table 4: Values of the scaling integrals dHJ)-© before rescaling and after rescaling, for m = 0. 



777 = 1 


B = 6 


B = 6 


B = 9 


B = 9 


B = 10 


B = 10 


B = 11 


B = 11 


II 


0.00882 


-8.15 10" 7 


-0.01094 


-6.73 10~ 9 


0.00975 


-2.61 10~ 9 


-0.00357 


6.74 10" 9 




0.00882 


-7.80 10- 7 


-0.01094 


-1.45 10" 9 


-0.01734 


7.01 10" 10 


-0.00357 


-2.06 10~ 8 


x 3 


-0.01765 


-2.27 10" 9 


0.02188 


-9.31 lO" 10 


0.00759 


4.16 10- 10 


0.00715 


-6.14 10" 9 



Table 5: Values of the scaling integrals (E])-® before rescaling and after rescaling , for m = 1. 
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- 2% differences in energy between the approximate Skyrmions before rescaling and the 
exact Skyrmion solutions. That implies that rescaling accounts for only a small part of the 
required change of the field needed to reach the exact solutions. Note also that nuclei with 
even baryon numbers have larger binding energies than those with odd ones. Encouragingly, 
Skyrmions with even baryon number generally have lower energy per baryon than Skyrmions 
with odd baryon number. This is seen in our results for B = 9, 10 and 11 presented in Table 

El 

Finally, baryon density isosurfaces are used to visualize the rescaled Skyrmion configu- 
rations and these are presented in Figure [TJ The shapes of the rescaled configurations are 
not very different from the unrescaled ones, since the values of the scaling parameters are 
close to 1. This can be seen in Figure El where plots of the unrescaled and rescaled baryon 
density isosurfaces of the 5 = 11 Skyrmion with massless pions are presented. In order to 
illustrate more clearly the deformation of the Skyrmions under rescaling we show in Figure 
El the B = 6 massless Skyrmion with the true rescaling parameters changed to the values 
Ai = A2 = 0.9 and A3 = I/A1A2, and the B = 9 massless Skyrmion with the rescaling 
parameters changed to Ai = A 2 = 1.1 and A 3 = l/AiA 2 . 




Figure 2: Isosurfaces of the unrescaled (left) and rescaled (right) baryon density for the 
B — 11 Skyrmion with m = 0. Each corresponds to a value of B = 0.035. 
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♦ # 



Figure 3: Isosurfaces of the true rescaled (left) and exaggerated rescaled (right) baryon density 
at a value of the baryon density which is one quarter of the maximum (i.e. 0.25 £> max ), for 
the B = 6 and the B = 9 Skyrmion with m = 0. 
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3 Conclusions 



We have started with a Skyrme field configuration of low but not minimal energy (JTJ), for se- 
lected values of B. This is constructed using the rational map ansatz. We have then rescaled 
the field, obtaining the energy (pjj), and have minimized this numerically, thereby finding a 
lower energy field configuration with the same symmetry. The rescaled field satisfies the 
scaling identities flS])-©, or equivalently (|9l-(fTTT). and is therefore closer to a true Skyrmion 
solution. In particular, we learn from the values of the rescaling parameters whether the 
true Skyrmion is more oblate or prolate than the rational map approximation suggests. This 
change of shape affects (slightly) the moments of inertia, which in turn will affect estimates 
of the energy spectrum of the quantized rotating Skyrmion, modelling real nuclei. 

The B = 6 configuration shrinks along the X3-axis and stretches in the {x\, X2)-plane. It 
has been observed previously that the rational map approximation to the B = 6 Skyrmion 
is slightly prolate, and rather too much so, in the sense that its classical electric quadrupole 
moment is positive, and rather too large to match the measured quadrupole moment of 
Lithium-6 [11] . The rescaling we have found makes the approximate Skyrmion less prolate, 
and this is better for fitting the quadrupole moment (although we have not recalculated it). 
Unfortunately this conclusion is not totally convincing, because if one compares the moments 
of inertia of the approximate B = 6 Skyrmion obtained using the rational map ansatz [TT] 
with the moments of inertia of the exact solution [9] , then it appears that the exact solution 
is more prolate. 

The B = 9 configuration stretches along the a^-axis and shrinks in the (x\, a^-plane, so 
it becomes more prolate. Similarly the B — 11 configuration stretches along the 23-axis and 
shrinks in the (xi, a^-plane, becoming prolate under rescaling. These are novel insights, and 
suggest in particular that the B = 11 solution is not closely related in shape to the known 
B = 12 solution with symmetry [9], which is significantly oblate. For B = 10 the three 
scaling parameters are independent, because the field configuration is triaxial, which agrees 
with the structure of the exact B = 10 Skyrmion (with m = 1) found numerically. 
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